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Spine-sheath layer radiative interplay in subparsec-scale jets and 
the TeV emission from M87 
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ABSTRACT 

Simple one-zone homogeneous synchrotron self-Compton models have severe difficulties is 
explaining the TeV emission observed in the radiogalaxy M87. Also the site the TeV emission 
region is uncertain: it could be the unresolved jet close to the nucleus, analogously to what 
proposed for blazars, or an active knot, called HST-1, tens of parsec away. We explore the 
possibility that the TeV emission of M87 is produced in the misaligned subpc scale jet. We 
base our modelling on a structured jet, with a fast spine suiTounded by a slower layer. In this 
context the main site responsible for the emission of the TeV radiation is the layer, while the 
(debeamed) spine accounts for the emission from the radio to the GeV band: therefore we 
expect a more complex correlation with the TeV component than that expected in one-zone 
scenarios, in which both components are produced by the same region. Observed from small 
angles, the spine would dominate the emission, with an overall Spectral Energy Distribution 
close to those of BL Lac objects with a synchrotron peak located at low energy (LBLs). 

Key words: galaxies: active-galaxies: jets-galaxies: individual: M87-radiation mechanisms: 
non-thermal. 



1 INTRODUCTION 

An increasing number of extragalactic objects is detected at ener- 
gies greater than 100 GeV. The list of published sources comprises 
20 objectiO: as expected (e.g., Costamante & Ghisellini 2002), the 
bulk of them (17) belongs to the class of Highly Peaked BL Lac 
objects. In fact, the position of the synchrotron peak, usually lo- 
cated in (or close to) the X-ray band assures the existence, in 
these sources, of electrons with extremely large Lorentz factors 
(7 = 10^-10^), a key ingredient for the emission in the TeV band 
via inverse Compton (IC) scattering. The remaining three sources 
are BL Lac itself, which belongs to the LBL class, (Albert et al. 
2007), 3C279, a FSRQs, (Teshima et al. 2007) and the nearby (16 
Mpc, Tonry 1991) radiogalaxy M87 (Aharonian et al. 2003, 2006). 

Already suggested as a possible source of high-energy radia- 
tion (e.g. Bai & Lee 2001) based on its similarity with BL Lac ob- 
jects (Tsvetanov et al. 1998), M87 was discovered as a TeV source 
by the HEGRA array of Cherenkov telescopes (Aharonian et al. 
2003) and it has been recently confirmed by H.E.S.S. (Aharonian 
et al. 2006). Due to the limited spatial resolution it is not possi- 
ble to identify the emission region of this radiation. However, the 
sensitivity of H.E.S.S. allowed the detection of variability on short 
timescales (~ days), suggesting a compact emission region. 

Three main sites have been proposed as emission regions for 
the TeV radiation: the resolved jet (in particular the so-called knot 
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HST-1); the unresolved base of the jet (in analogy with blazars); 
the vicinity of the supermassive (M = 3 x lO'' Mq, Marconi et al. 
1997) black hole (BH). 

HST-1 is a quite interesting jet feature located at 60 pc (pro- 
jected) from the core of of M87. Superluminal motions with ap- 
parent speeds up to ~ 6c were observed with HST (Biretta et al. 
1999). Continuous multifrequency coverage has shown spectacular 
activity from this knot in the past years (Perlman et al. 2003; Harris 
et al. 2006, 2007; Cheung et al. 2007). In particular, the X-ray flux 
increases by a factor of ~ 30 from 2000 to 2005 and since then it 
is steadily decreasing. These variations are accompanied by similar 
changes in the optical and in the radio bands and by the appear- 
ance of new superluminal components. This extreme phenomenol- 
ogy is described by Stawarz et al. (2006) on the assumption that 
HST-1 marks the recoUimation shock of the jet. As such, HST-1 is 
thought to be a rather efficient particle accelerator and thus a possi- 
ble source of intense TeV radiation. However, the small variability 
timescales observed by H.E.S.S. put strong limits to the size of the 
source(ii < cAtS ~ 5x 10^^(5 cm, where.? = [r{l- f3 cosO)]-^ 
is the relativistic Doppler factor, 6 is the viewing angle and ~ 2 
days), which seems difficult to accomplish in this scenario (e.g. 
Neronov & Aharonian 2007; Levinson 2007). 

The rapid variability could be easily reproduced if the TeV 
emission is produced close to the base of the jet, in the region as- 
sociated to the emission from blazars. However, standard one-zone 
leptonic models face difficulties in reproducing the observed spec- 
tral energy distribution (SED) of the core including the TeV data. A 
possible solution advocates that the emission comes from a decel- 
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erating jet (Georganopoulos et al. 2005) . Alternatively, the emis- 
sion can be produced by high-energy protons (Reimer et al. 2004). 
Going further in towards the BH, TeV photons could be emitted 
through IC by relativistic pairs produced by electromagnetic cas- 
cades in the BH magnetosphere (Neronov & Aharonian 2007). 

The determination of the site producing the high-energy emis- 
sion is rather important. Indeed, if the emission site will be even- 
tually identified as knot HST-1, this would have a broad impact on 
the current view of relativistic jets (see the discussion in Cheung 
et al. 2007). It is thus extremely important to investigate whether 
the standard view (or minimal variations of it) for the inner jet is 
able reproduce the observed phenomenology. If not, it would be 
mandatory to explore the other alternatives. 

As discussed in Ghisellini, Tavecchio & Chiaberge (2005, 
hereafter Paper I), there is compelling evidence supporting the view 
that jets of BL Lac objects are structured at blazar scales (at dis- 
tances ~ 10^^ cm from the BH), with a fast core (spine) surrounded 
by a slower sheet {layer). In Paper I we showed that the radia- 
tive coupling between the spine and the layer can naturally account 
for the deceleration inferred for the jet between the blazar and the 
VLBI scale. Furthermore, one component sees the radiation of the 
other boosted, since the relative velocity can be relativistic. This 
enhances the IC emission of both the spine and the layer 

A direct consequence of such a structure is that while at small 
viewing angle (as in the case of blazars) the emission is dominated 
by the boosted spine emission, at large observing angles (6* > 45°, 
typical for radio-galaxies) the emission from the spine would be 
suppressed, while the layer, characterised by a broader beaming 
cone, could substantially contribute to (and sometimes dominate) 
the overall emission. At intermediate angles both components can 
significantly contribute. Following this line, in this letter we explore 
an alternative interpretation for the TeV emission of M87, assuming 
that the TeV radiation is produced in the layer of the misaligned jet. 
We use Ho = 70 km Mpc"\ = 0.7, = 0.3. 



2 THE SPINE-LAYER SCENARIO 

2.1 Difficulties witli tlie one-zone models 

The SED of the core of M87 is reported in Fig.[T](grey open points, 
see Paper I for references) together with the H.E.S.S. spectra for 
2004 (open blue squares) and 2005 ( open red triangles) as reported 
by Aharonian et al. (2006). The X-ray spectrum (as observed in 
July 2000) is taken from Balmavarde, Capetti & Grandi (2006). 
Note that the slope [a — 1.3 ± 0.1; F{i/) cx f'"] is intermediate 
between that reported in Marshall et al. (2002) (a = 1.47 ± 0.08, 
see also Donate, Sambruna & Gliozzi 2004) and that of Wilson & 
Yang (2002) (a — 1.17 ± 0.1). The X-ray spectral shape was not 
determined in 2004 and 2005, and for simplicity we assume that it 
is the same of 2000, with a normalisation that is a factor 3 (2004, 
blue) and 4 (2005, red) greater. This reproduces the increase of the 
X-ray flux of the core recorded through the continuous monitoring 
of Chandra (Harris et al. 2007). 

The data describe a bump peaking at ^ 10^* Hz and extend- 
ing into the X-ray band, while the TeV emission clearly belongs to 
a second component. The SED closely resembles that usually ob- 
served in blazars, in which the first component is synchrotron emis- 
sion from relativistic electrons and the second peak is produced 
through IC emission (e.g., Ghisellini et al. 1998). For M87 there is 
no evidence of powerful sources of soft photons in the nucleus and 
hence the high-energy emission should be dominated by the Syn- 
chrotron self-Compton emission (SSC; e.g., Maraschi et al. 1992). 



However, a simple one-zone synchrotron-SSC model cannot eas- 
ily describe the observed SED. The reason is that in order to have 
the synchrotron component peaking in the IR region (fs — 10^* 
Hz) and the SSC one close to the TeV band (uc — 10^* Hz, as 
required by the relatively flat TeV spectra) one is forced to assume 
a rather low magnetic field and an unreasonably large Doppler fac- 
tor Following Tavecchio, Maraschi & Ghisellini (1998) and defin- 
ing the synchrotron and SSC peak luminosities, Ls = h'sL{vs) and 
Lc = vcL{vc) and the variability timescale At, one obtains: 

5 ^ 500 L\%,,L-%^ ujl, u'Jl, M-'J^ (1) 

Here Q — 10"^ Q^,, in cgs units, except for At, that is in days. The 
adopted numerical values are appropriate for M87. 

This is a general problem faced by any synchrotron-SSC 
model using a single emission region. A direct way to overcome 
this problem is to "decouple" synchrotron and the IC components, 
assuming that the two components are produced by two different 
regions. Such a "two-zones" model would allow us to reproduce 
the low-frequency bump (emitted by electrons with relatively low 
energy) with one emission component and the TeV emission with 
a separate source with high-energy electrons. This opportunity is 
directly offered by the spine-layer scenario. In this case the de- 
beamed spine (or, alternatively, the layer) could be responsible for 
the low-energy bump and the layer (or the spine) of the TeV com- 
ponent. 



2.2 Tlie model 

We focus on the case of a sub-pc scale emission region, due to 
the observed inter-day variability of the TeV emission, difficult to 
explain if HST-1 is responsible for this emission. 

We are aware that we are almost doubling the number of the 
free parameters (one set for the spine, one set for the layer), and 
therefore we pay the increased flexibility with an increased uncer- 
tainty concerning the value of the parameters themselves and a loss 
of predictive power. On the other hand, our main aim at this stage is 
not much to find whether/what specific physical conditions (value 
of the magnetic field, beaming factor, number of particles, and so 
on) are univocally determined by the model, but to show that one 
can attribute the radiation we see to the the sub-pc jet (with "minor" 
modification of the simplest model) without radically changing our 
ideas about the emission models of blazars. Having said that, we 
will nevertheless do our best to find the most reasonable set of pa- 
rameters. This is why we will not only try to fit the SED of M87, 
but we will worry about what kind of blazar would M87 be if seen 
at small angles. If the resulting SED would be too odd with respect 
to the known blazar SED, we will discard the found parameter set. 

We refer to Paper I for a complete description of the model. 
Briefly, we assume that the jet comprises two regions: i) the spine, 
assumed to be a cylinder of radius R, height Hs (as measured in 
the spine frame) and in motion with bulk Lorentz factor F^; ii) the 
layer, modelled as an hollow cylinder with internal radius ii, ex- 
ternal radius i?2, height Hi (as measured in the frame of the layer) 
and bulk Lorentz factor F;. Each region is characterised by tangled 
magnetic field with intensity Bs, Bi and it is filled by relativistic 
electrons assumed to follow a smoothed broken power-law distri- 
bution extending from 7min to "jcut and with indices ni, 712 below 
and above the break at 76 . The normalisation of this distribution 
is calculated assuming that the system produces an assumed syn- 
chrotron luminosity Lsyn (as measured in the local frame), which 
is an input parameter of the model. As in Paper I we assume that 
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Figure 1. SED of the core of M87 (grey open points) together with the 
H.E.S.S. spectra taJcen in 2004 (open blue squai'es) and 2005 (open red tri- 
angles), from Ahai'onian et al. (2006). The grey bow-tie reports the X-ray 
spectrum as measured by Chandra in 2000 (from Balmaverde et al. 2006). 
We reproduced the increased X-ray emission of 2004 and 2005 assuming 
the same slope and a larger normalisation (blue and red bow-ties). The lines 
report the emission from the spine and from the layer for the two states as 
calculated with the parameters reported in Tab. 1 . For comparison the green 
line indicates the sensitivity of GLAST (5cr, 1 year, converted in luminosity 
assuming the distance of M87). The upper black line reports the emission 
from the spine that would detect an observer located at an angle of 6° from 
the jet axis (the corresponding emission from the layer is well below that of 
the spine and is not reported for semplicity). For comparison, we show the 
data of BL Lac (filled orange symbols, from Ravasio et al. 2002, TeV data 
from Albert et al. 2007). See text for details. 

Hi > Hs . The seed photons for the IC scattering are not only those 
produced locally in the spine (layer), but we also consider the pho- 
tons produced in the layer (spine). Hereafter we will indicate this 
component as External Compton (EC). 

To properly calculate the observed emission, it is necessary to 
take into account the different beaming patterns associated to the 
two emission components. This is because, in the comoving frame 
of the spine (layer), the photons produced by the layer (spine) are 
not isotropic, but are seen aberrated. Most of them are coming from 
almost a single direction (opposite to the relative velocity vector). 
In this case the produced EC is anisotropic even in the comoving 
frame, with more power emitted along the opposite direction of the 
incoming photons (i.e. for head-on scatterings). We therefore must 
take into account this anisotropy in the comoving frame when trans- 
forming in the observer frame. Consider also that this applies only 
to the EC radiation, while the synchrotron and the SSC emissions 
are isotropic in the comoving frame. The EC radiation pattern of the 
spine will be more concentrated along the jet axis (in the forward 
direction) with respect to its synchrotron and SSC emission. On the 
contrary, the EC radiation pattern of the layer will be less concen- 
trated along the forward direction (with respect to its synchrotron 
and SSC) because, in the comoving frame, the layer emits more EC 
power in the backward direction (e.g. back to the black hole). The 
corresponding transformations are fully described in Paper I, and 



we briefly summarise them in the following. We here follow Der- 
mer (1995), assuming that, due to the strong aberration, in the rest 
frame of the spine (layer) all the photons, assumed to be isotropic in 
the layer (spine) frame, can be considered as coming from a single 
direction, i.e. along (opposite to) the jet axis. This approximation 
greatly simplifies the calculations, allowing us not to consider the 
detail of the geometry in the transformations. 

The beaming of the synchrotron and SSC emission for both 
regions follows the usual S'l'^" (spine) and (5f^° (layer) pattern, 
while the patterns of the EC components are described by more 
complex shapes. In particular, the EC from the spine will have 
a rather narrow pattern given by (J^^^^^f^", where 5s,i is the 
Doppler factor of the spine as measured in the layer frame. The 
first term {5'^^'^°') describes the same pattern as the EC radiation 
produced when scattering the seed photons produced by the broad 
line regions of powerful blazars, as firstly pointed out by Dermer 
(1995). It can be understood by going in the reference frame where 
the seed photons are isotropic (in this case this means comoving 
with the layer). In this frame we see the EC from the spine boosted 
by the "Dermer" pattern 5^^'^°' ■ Finally, we go to the observer 
frame using the standard 5f term. Analogously, the amplifica- 
tion of the EC from the layer is oc (5f+^°(5f+°'. The net result is 
that the EC emission from the layer is less boosted in the observer 
frame than the corresponding EC emission from the spine. 

Fig. [2] illustrates this point, showing the amplification factors 
of the different components as a function of the viewing angle 
(from 0° to 45°), calculated assuming Vs = 12 and F; = 4 (the 
values we used for M87, see below). Solid lines refer to the spine, 
dashed lines to the layer. Black and red lines correspond to the syn- 
chrotron/SSC and EC amplification factors, respectively. For the 
spine, the EC amplification is larger than the Synchrotron/SSC one 
for small angles, while at large angles the latter dominates. In this 
sense one can say that the EC emission has a narrower pattern than 
the Synchrotron/SSC emission. The opposite occurs for the layer: 
at small angles the synchrotron/SSC amplification is larger then 
the EC one, which instead becomes more important at large an- 
gles. Note, moreover, that the EC/layer emission has its maximum 
at ~10°, while all the other curves peak at 0°. To understand this, 
consider the emission pattern in the frame comoving with the layer. 
In this frame the seed photons produced by the spine are coming 
along the black hole-layer direction: as a consequence the relativis- 
tic electrons will produce more EC radiation in the opposite (layer- 
black hole) direction. This is the radiation that will be boosted the 
least when we transform the intensity in the observer frame (and 
the opposite for radiation along the outward radial direction: in the 
comoving frame it is the weakest, but it is boosted the most when 
going to the observer frame). This "compensation" between the in- 
trinsic anisotropy as seen in the comoving frame and the beaming 
effect means that the maximum flux will be observed not at 0° , but 
at a larger angle. At the extreme, if the layer is not moving, the 
maximum will be observed at 180°. This characteristic is partic- 
ularly relevant for sources seen at large angles: for those, the EC 
produced by the layer becomes relatively more important. 



2.3 Results 

We have two possibilities to reproduce the SED with this model: 
i) either we assume that the low energy bump is produced by the 
layer and the TeV emission originates in the spine or ii) the layer 
produces the TeV radiation and the spine mainly contributes at low 
frequencies. The first choice is disfavoured in the present case. In 
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Figure 2. Amplification factors for tlie different emission components as a 
function of the viewing angle. Solid (dashed) lines refer to the spine (layer), 
while black (red) lines show the to synchrotron/SSC (EC) factors. We as- 
sume Fa = 12 and F; = 4, the values we used for the case of M87, and 
a = 1. The two vertical lines indicate the two angles assumed for the cal- 
culation of the SEDs, 6° and 18° . See text for discussion. 



fact, to produce TeV photons, the typical Lorentz factor of the elec- 
trons in the spine must be large (7 ~ 10'') and, assuming typical 
magnetic fields of order 0.1-1 G, the corresponding synchrotron 
emission would peak in the hard X~ray band. Under this conditions 
the SSC emission would be deeply in the Klein-Nishina regime 
and thus strongly depressed and the TeV emission would be dom- 
inated by the EC component. If the spine SED is dominated by 
the EC component at large (18°) angles, it will be even more so at 
small angles (see Fig. 2), with LtcV / ix 3> 20. Instead, in TeV 
BL Lacs, the synchrotron and the IC components show compara- 
ble luminosities [with few exception, most notably lES 1426-1-428 
(Aharonian et al. 2002), although the TeV luminosity depends on 
the still unclear amount of intergalactic absorption]. On the other 
hand, as we show below, the choice ii) leads to a predicted SED 
for the beamed counterpart of M87 closely resembling those ob- 
served in LBL sources and thus we decided to adopt this view in 
the modelling. 

The SEDs obtained for the core of M87 with the model are 
shown in Fig. [T] for both 2004 and 2005 and the parameters are 
reported in Table 1 . 

The model shown is intended as an example of the possible 
parameters giving a reasonable shape for the SED. A whole fam- 
ily of acceptable solutions, each one with different values for the 
parameters, is allowed. However, the joint request to reproduce the 
data and to have a spine SED resembling that of a BL Lac at small 
angles provides some constraints. Since the layer has to produce 
TeV radiation, the Lorentz factor of the electrons at the peak is con- 
strained to be of the order of 7;, ~ lO'', independently on the com- 
ponent (SSC or EC) dominating the high-energy emission. This, in 
turn, fixes the peak of the synchrotron component close to or within 
the X-ray band for the chosen magnetic field. Smaller B-values 
are allowed, but not larger, since in this case the layer synchrotron 
emission would overproduce the observed data. For the spine, as 
already discussed, most of the constraints come from the request 
that the Lorentz factor, the size and the magnetic field are similar 
to those typically inferred for (LBL) blazars. 

The model slightly overproduces the emission in the optical 
band. However, since the optical and the X-ray emission are clearly 



correlated (Perlman et al. 2003) we expect that high X-ray fluxes 
correspond to high optical states. Note also that we predict quite 
bright GeV emission in this period, exceeding the EGRET upper 
limit by a factor of 3 and easily detectable with GLAST. The high- 
energy component of the layer is dominated by the EC emission 
from the beamed (in the layer frame) IR-optical photons from the 
spine. On the other hand, since the synchrotron component of the 
layer peaks in the X-rays, the EC emission of the spine is strongly 
suppressed by the effect of the Klein-Nishina cross section and 
hence the high-energy component is mainly SSC. The variability 
of the TeV flux from 2004 to 2005 is reproduced by increasing 
the luminosity radiated by the layer (and thus the corresponding 
density of the electrons) and slightly decreasing the magnetic field. 
Note also that the small size of the emitting region allows the rapid 
changes such as those detected by H.E.S.S. (At ~ 2 days). 

A supplementary contribution to the IC emission from both 
the spine and the layer could come from the scattering of the ex- 
ternal IR radiation field, possibly associated to the putative low- 
efficiency accretion flow. The luminosity of the core of M87 in the 
infrared band is L « 10*^ erg s~^ (Perlman et al. 2001). This 
is an upper limit to the luminosity of the accretion flow, since the 
observed power law spectrum supports the idea that it is mainly as- 
sociated to the emission from the jet (Perlman et al. 2001). We can 
set an upper limit to the importance of this component by assum- 
ing L = 10'*^ erg s~^ for the accretion luminosity, produced in a 
region whose size is of the order the distance of the active region 
of the jet (around ~ 10^^ cm). This maximizes the importance of 
this external component: for smaller radii the accretion IR radiation 
would be strongly deboosted in the jet frame, while for larger radii 
it will be more diluted. Even so, this external component has a radi- 
ation energy density smaller (by more than an order of magnitude) 
of the one already present in the spine and in the layer. 

Inspection of Fig.[T]reveals that we nicely reproduce the slope 
of the 2004 TeV spectrum, while the 2005 spectrum appears to be 
harder than the model (although the relatively large error bars pre- 
vent a firm conclusion). The steep slope of the model continuum in 
this band is the result of the 77 absorption of TeV photons interact- 
ing with optical-IR photons of the spine and thus cannot be simply 
avoided. 

In Fig. [T] we also report (black line) the emission from the 
spine relative to 2004 as would appear to an observer almost 
aligned {6 — 6°) with the jet. For comparison we report (orange, 
filled symbols) the multifrequency SED for BL Lac itself (the pro- 
totype of BL Lac objects), collected during different campaigns 
(see Ravasio et al. 2002 for references) and reporting also the TeV 
spectrum recently detected by MAGIC (Albert et al. 2007). 



3 DISCUSSION 

We propose that the TeV emission detected from M87 originates 
in the structured jet at blazar scale. The most likely possibility is 
that the emission from the spine reproduces the low energy compo- 
nent, from radio to X-rays, while the layer contributes to the bulk 
of the TeV radiation. With this choice, the beamed counterpart of 
M87 observed at small angle would have a SED closely resembling 
those of Low-Peaked BL Lac objects. Correspondingly the spine is 
characterised by physical parameters close to those usually inferred 
for those sources (75 ~ 10^, _B ~ 1 G). The layer, instead, would 
be characterised by a low magnetic field (B ~ 0.1 G) and large 
peak Lorentz factors (7;, ~ 10®). It is possible that particles in the 
layer are energized by turbulent acceleration (e.g. Stawarz & Os- 
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Table 1. Input parameters of the models for the layer and the spine shown in Fig.[T] All quantities (except the bulk Lorentz factors F and the viewing angle 8) 
are measured in the rest frame of the emitting plasma. The external radius of the layer is fixed to the value R2 = 1.2 X R. The corresponding values of the 
electron density are: 108.3 cm""^ and 140.8 cm~^ (spine low and high); 0.74 cm~^ and 2.04 cm~^ (layer low and high). 



trowski 2002). In this case the expected distribution would have a 
pile-up at the energy where acceleration and loss processes are in 
equilibrium (e.g. Katarzynski et al. 2006). Suggestively, in our con- 
ditions (losses dominated by IC) the expected equilibrium Lorentz 
factor would be close to 7 ~ 10^ for B ~ 0.1 G. 

Note that our conclusions are rather different from that of 
Georganopoulos et al. (2005) who discussed a similar model re- 
lying on the possibility that the jet experiences a strong decelera- 
tion at the blazar scale. In their model the high-energy emission is 
produced in the innermost and faster portion of the jet, while the 
slow external portions mainly contributes at low frequency. Under 
these conditions the beamed counterpart of M87 would be a TeV 
BL Lac. However, due to the different beaming patterns or the syn- 
chrotron and IC components (analogously to our case), the result- 
ing spectrum would be characterised by a strong dominance of the 
TeV component with respect to the X-ray one, contrary to what is 
usually observed in the known TeV BL Lacs. 

In our model the optical and the X-rays are produced by a 
region different from that responsible for the TeV emission. There- 
fore a strict correlation between low energy and TeV emission is 
not directly required (even if it is possible). Analogously, the MeV- 
GeV emission should originate in the spine. Therefore we expect 
that the emission possibly detected by GLAST will not exactly fol- 
low the TeV component. 

A potential weak point of our model concerns the slope of the 
TeV spectrum. The slope is mainly dictated by the absorption of 
TeV photons by the dense optical radiation field and, in this sense, 
it is rather robust and does not strongly depend on the underlying 
slope of the intrinsic TeV spectrum. In particular, hard spectra such 
as that recorded in 2005 are quite difficult to achieve in this scheme. 
A possibility to avoid important absorption of gamma-rays would 
be to enlarge the emission region, hence diluting the target photon 
density and decreasing the absorption optical depth. However the 
increase of the source size is limited by the observed short vari- 
ability timescales. Further observations with increased sensitivity 
could provide stringent constraints to our interpretation. 
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